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I. INTRODUCTION 



A. EARTHQUAKE PREDICTION THEORY 

The desirability of being able to predict when and where 
an earthquake will occur becomes immediately apparent in 
light of the recent devastation that occurred in Mexico 
City. The potential for death if a catastrophic earthquake 

should occur along the San Andreas fault is estimated to be 
between 3,000 and 13,000 lost lives. The lower figure 
applies for a quake that hits at 0230 and the higher figure 

applies for a quake that occurs at 1630 on a weekday. These 

predictions are based on the facts that wooden, single 
family homes are much more resistant to earthquake struc- 
tural damage than are many of the older, unreinforced 
masonry buildings. In May, 1 983, 30 of the 40 unreinforced 

masonry structures in Coalinga, Ca. , either partially or 
totally collapsed. In Los Angeles, 8,000 similar buildings 

exist. The probability of a magnitude 8.3 earthquake occur- 

ring along the southern San Andreas fault within the next 30 
years is estimated to be 50 percent. [Ref. 1 : p . 35] 

The primary cause of earthquakes is the movement of the 
litnosphere, the outer layer of the Earth consisting of 
twelve rigid plates riding on the supporting mantle. This 
movement rearranges continents and forms and destroys both 
mountains and oceans. The driving energy for this poten- 
tially large scale destructive mechanism is the convective 
circulation of the mantle. The mantle consists of hot, 
solid rock. In time measured on a geological scale, the 
mantle is sufficiently hot that it can flow and be easily 
deformed. Hot magma from the mantle wells up between sepa- 
rating plates at midocean ridges forming new lithosphere. 
The movement of the cold, brittle plates is primarily along 
Boundary faults and is controlled by friction. Earthquakes 



slip and elastically 



occur when these plates suddenly 
rebound. [Ref. 2] 

Figures 1.1 and 1.2 [Ref. 1:p. 38] show the typical 

cycles of earthquakes as explained by the theory of plate 
tectonics. In figure 1.1 the hot mantle 10-15 km deep is 
plastically deforming and stress is accumulating at asperi- 
ties, or lock points, in the cooler crust. At the moment of 
rupture an asperity breaks and crustal rupture spreads 
causing earthquakes. The slippage is halted when it encoun- 
ters another asperity in approximately 100 seconds. 
Aftershocks may occur for the next several months. 
Additional slippage at or near the surface may or may not 
occur. In Figure 1. 2 a second major earthquake has occurred 
50 or 100 years after the first. The cycle begins again 
after about 200 years when the deep crust displacement 
imposes new stress on the upper brittle crust. 

To date, earthquake prediction has been on a long-term 
basis normally expressed as a probability of quake occur- 
rence during the coming scores of years. What has been 
lacking is a reliable, short-term prediction process. 

Central to any short-term technique is geologic phenomena 
during the period immediately preceding the main earthquake. 
Before the arrival of a major shock, swarms of many small 
earthquakes are often noted or sometimes accelerated defor- 
mation of the crust is observed. However, neither one of 
tnese two predictors is reliably present in all cases. 
[Ref. 1 : p. 38] 

Systematically collected data on the generation of elec- 
tromagnetic energy prior to, during, or after an earthquake 
currently does not exist. Evidence of radio frequency (RF) 
emissions as a precursor to quake activity appears impul- 
sively suggesting these waves are generated over a broad 
spectrum of frequencies [Ref. 3]. Citizen band and amateur 
radio operators in the Hollister, Ca. , area reported an 
"increase in background noise level preceding earthquake 




Figure 1.1 First Portion of an Earthquake Cycle 




Figure 1.2 Second Portion of an Earthquake Cycle 



activity" at 27 MHz and at the six and two meter bands. 
This increase in noise was on the order of 10 dB and 
preceded a quake by 12 to 24 hours. At one hour to 15 
minutes before a quake, noise levels returned to normal. 
The appearance of earthquake lights [Ref. 4] and reception 
of 18 MHz radio noise [Ref. 5] are further examples of a 
possible connection between RF energy and earthquakes. 



Various theories explaining the actual mechanism of RF 
generation at the fault boundary include impact of 
dissimilar metals, shock induced polarization and fracture 
of guartz-bearing rock [Ref. 6]. Because of the strong 
implication that earthquakes generate electromagnetic 
energy, this thesis will present a satellite design for a 
space-based system as a means of efficient quake detection 
and rapid data dissemination. 

B. SEACE-BASED EARTHQUAKE DETECTION 

The exploitation of outer space as a tool in earthquake 
prediction based on electromagnetic detection can be divided 
into two major scenarios. The first is installing ground- 
based sensors that detect the quake generated RF energy and 
then relay data to an orbiting satellite. The satellite 
would function as a central collection point and could 
incorporate features such as data processing and earthquake 
pattern recognition with automatic alarm transmission to key 
Earth stations. The number of ground-based sensors to 

provide adequate coverage can not be reliably estimated 
because the strength of the RF energy is not known but is 
being researched as a portion of this thesis. As the 

strength increases, the distance between sensors can 
increase and their number decrease. 

The second scenario is direct detection of earthquake 
generated RF energy by an orbiting satellite. Ail the 

previously mentioned features could be included. Direct 
satellite detection would require a sufficiently strong RF 
signal to overcome the power spreading loss incurred from 
the large distance traveled. The major advantages over the 
first scenario are worldwide coverage with elimination of 
the need for ground-based RF sensors. On an international 
level, frequencies could be allocated to quake detection to 
minimize possible man-made interference with the RF signal. 
Because ox the advantages, this second scenario will be the 
basis for the proposed satellite design. 
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C. FCCUS OF STUDY 

A proposed design for a space based sensor to predict 
and detect earthquakes is presented- A free standing BF 
noise measurement and recording system is described to 
research the correlation between earthquakes and increased 
background electromagnetic noise at 30.45 MHz and 150.75 
MHz. This thesis continues the research, construction and 
installation of equipment described in Reference 6. 



II. BACKGROUND 

A. ORBITAL PARAMETERS 

1 . Space Coordinate System 

In order to understand the background information 
that will be used in subsequent portions of this thesis, it 
is necessary to first establish the coordinate system that 
describes a satellite's position in space. An ephemeris is 
a tabulation of the predicted or observed positions occupied 
by the satellite. 

The coordinate system to be used is shown in Figure 
2.1 [Ref. 7]. The large, flat rectangular plane is an 
extension of the plane of the equator and is called the 
equatorial plane. Superimposed on the equatorial plane is 
the Cartesian coordinate system indicated by I, J and K with 
the origin at the Earth's center. These coordinates do not 
rotate but are fixed in space with the I axis pointing in 
the direction of the vernal equinox. 

The satellite's motion about the Earth is on the 
circumf erence of the orbital plane. An elliptical orbit is 
shown but a satellite can also have a circular orbit around 
the Earth. The angular momentum or spin vector (h) for tae 
satellite orbit is determined by the right-hand rule for a 
south to north eguator crossing and is perpendicular to the 
orbital plane. Inclination (i) is measured in degrees from 
K to the spin vector h. By definition, inclination can vary 
from a minimum of zero degrees to a maximum of 180 degrees. 

The eccentricity vector (e) starts from the origin 
and points in the direction of perigee or the periapsis 
direction. Perigee is the satellite's closest point of 
approach to the Earth. The line of nodes (n) is the 
intersection of the equatorial plane and the satellite's 
orbital plane. The angle between the line of nodes and the 
eccentricity vector is called the arguement of perigee and 
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Figure 2.1 Classical Orbital Elements 
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is measured in the direction of the satellite motion. 
Right ascension of the line of nodes (n) is the angle 
measured in the equatorial plane from the vernal equinox 
direction to the line of nodes "n. The epoch angle (v 0 ) is 
the angle measured from e - to the vector from the center of 
the Earth to the satellite (r„ ). The range of the epoch 
angle is from zero degrees to 360 degrees. 

2. Equations f o r an Ellipse 

Many satellites orbit the Earth in an elliptical 
manner or use portions of an elliptical orbit during an 
orbital transfer from one circular orbit to another circular 
orbit of a different radius. Because the equations of an 
ellipse are an integral portion of satellite orbitology, the 
pertinent relationships will be presented. Figure 2.2 shows 
the classical form of an ellipse with the small circle 
representing the Earth at one focus. 

The semi-major axis (a) is the horizontal length 
measured from the center of the ellipse to the outer edge. 
The semi-minor axis (b) is the corresponding vertical 
length. The eccentricity (e) is a measure of departure from 
a cirlce and is given by Equation 2.1 

e = SQET (1- (b/a) 2) (2. 1) 

Eccentricity can vary from a value of zero indicating b = a 

for a circle to a value of slightly less than one if b 

approaches zero. 

The radius of apogee (R a ) is the distance from the 
center of the Earth to the point of farthest distance from 
the Earth. The radius of perigee (R p ) is measured from the 
center of the Earth to the satellite's closest point of 

approach to the Earth. The relationship between R a and e 

can be expressed as 

R a = a (1 +e) (2.2) 
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Figure 2.2 An Ellipse Around the Earth 
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A similar formula Holds for R p and e. 

Rp = a(1-e) (2.3) 

B. ORBIT EQDATIOHS 

1 - Ci rcular and Elli pt ica 1 Orbit Vel ociti es 

Two celestial bodies, one in orbit around the other, 
comprise a two-body system and must obey Newton's law of 



universal gravitation 


given by 


F = GMm 
r 2 


(2.4) 



Referring to Figure 2.3, F is the force of attraction due to 
gravity between masses a and m separated by a distance r. 
The constant G is the universal gravitational constant and 
its value is 6.67 X 10~ 11 m 3 /k g sec 2 . For any body of mass 
m in a circular orbit with velocity v and at a distance r. 



the centripital force 


is given by 


F = m v 2 
r 


(2.5) 


The centripital force 


must equal the gravitational force or 



the orbiting body will not be in an equilibrium orbit. 

GJm = mv 2 (2.6) 

r 2 r 



Cancelling the common 
for v yields 


terms of m, one of the r's and solving 


v = 5 QRT (GM/r) 


(2.7) 



This v is the linear velocity of mass m in a circular orbit 
about mass M at a radial distance r. The time required for 
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m to complete one circular revolution (T c ) is found by 
dividing the distance traveled {2z-rr) by the velocity. This 
gives 



T c = 2r 7r (2.3) 

v 

= 2r 7T (2.9) 

SQRT (GM/C) 

= 2 tt r 3 / 2 (2. 10) 

SQRT (GM) 

Kepler's third law of planetary motion relates the 
motion of one planet to another. This relationship is inde- 
pendent of the shape of their orbits and requires that the 
squares of the periods or the time to complete one revolu- 
tion vary as the cubes of the semi-major axes. Then for m 

in an elliptical orbit about fl with a semi-major axis of a, 
the time for one revolution (T e ) is given by 

T e = 2 ttSl 3/2 (2. 1 1) 

SQRT (GM) 

The elliptical orbit velocity (V e ) is given by 

V e = SQRT[ Gil (2/r- 1 /a) ] (2.12) 

2 . Orbi tal T ran sfe r s 

It is often desired to change the omit of a satel- 
lite to place it either in a more advantageous relative 
position to other satellites or to alter its relative posi- 
tion to the Eartn and allow it better access no Earth gener- 
ated data. To perform an orbit modification the satellite 
must expend energy, normally in the form of burning fuel 
from an on-board rocket or releasing compressed gas tnrough 




Figure 2.3 Gravitational Force of Attraction 

thrusters. Rockets are used when a major change is desired 
because a large force can be quickly developed while 
compressed hydrazine gas powered thrusters are normally used 
for station keeping or very minor orbital adjustments. 

The two major factors that influence exactly how an 
orbit modification will be accomplished are fuel and time. 
If it is desired to achieve the new orbit as rapidly as 
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possible, then a large amount of fuel will be required to 
provide the necessarily large force and acceleration- Such 
a time constraint may exist in the case of rescue missions 
or anti-satellite weapons. If time is not the deciding 
criterion then slower, more fuel efficient transfers are 
possible. This has the major advantage of allowing more of 
the satellite's weight to be in instrumentation since less 
weight need be allocated to fuel. A fuel efficient transfer 
will be used for this thesis to maximize operational 
payload. 

Figure 2.4 depicts a satellite at position one in a 
counter-clockwise circular orbit about the Earth. It is 
desired to discharge a payload from the satellite so that 
upon completion of the maneuver, the payload will be in 
position two and the satellite will return to position one. 
The positions are separated by the angle u measured in 
radians. This same type of maneuver could be used if posi- 
tion one were to maneuver in order to rendezvous with a 
second satellite already at position two. To accomplish 
either task, the satellite at position one will accelerate 
into the elliptical omit shown. While position two 
continues on its circular orbit and may pass the original 
position one several times, the accelerated satellite will 
travel the elliptical path and rendezvous will occur when 
both satellite one and two meet at position one. Note that 

tne satellite at position one can not merely "speed up" to 
catch position two because by Equation 2.7 they would be at 
different distances from the Earth. It is necessary to know 
how much position one's velocity must be increased to place 
it into the correct elliptical orbit so rendezvous will 
occur . 



Since the satellites 
the time for one elliptical or 
for an integer number of circu 
circular orbit to account 



are to meet at position one, 
bit must be equal to the tiae 
lar orbits plus a portion of a 
for the initial angular 
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Figure 2-4 Circular Orbit Modification 
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difference. The integer number of orbits (K) can be 
selected to again either minimize time or fuel expended with 
more circular orbits indicating less fuel required. That is 

KT C + T c (1-u/2ir) = T e (2.13) 

Substituting Equations 2.10 and 2.11 into Equation 2.13 
giv es 

K2ttR c 3/2 + 27tR c 3/2 (1-u/2tt) = 2^3/2 (2.14) 

SQBT(GM) SQRT(GM) SQRT (GM) 

The circular orbit radius is now being represented by R c 
rather than just r to assist in equation development. The 
common term of 2 tt/SQ RT (GM) can be cancelled which leaves 

KR C 3/2 + R c 3/2 ( 1-U/27T) = a3/2 (2.15) 

Rearranging and solving for a gives 

a = R c [K + l-u/2 :r]2 /3 (2. 1o) 

Equation 2.16 gives an expression for a which is the 
semi-major axis of the required elliptical orbit. An 
expression is needed that relates the eccentricity of the 
transfer ellipse to K. This can be done by using Equation 
2.3 and realizing that at perigee R c = R p . Maxing the 
substitution yields 

a = a (1-e) [ K+l-u/2 fr]2/3 (2.17) 

Cancelling the a and rearranging gives 

1 = [ K + l-u/2 r]2 /3 - e[ K+l-u/2 tt]2/3 (2.18) 

Solving for e gives Equation 2.19 
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e = [ K+1 -u/2tt] 2/3 
[ K + l-u/2 tt] 2 



1 



( 2 . 19 ) 



Equation 2.19 now relates the necessary eccentricity 
of the transfer ellipse to both the number of circular 
orbits of position two and to their initial angular separa- 
tion. Equation 2.7 gives the circular orbit velocity and 
Equation 2.12 describes the elliptical orbit velocity- The 
velocity of the elliptical orbit at R p must be known to find 
the difference between the two velocities. Evaluating 
Equation 2.12 at r = R p gives 

V p = SQET[GM(2/R p - 1/a) ] (2.20) 

V p is the velocity at perigee for satellite one. Solving 
Equation 2.3 for a gives 

a = R p (2.21) 

1-e 

Substituting Eguation 2.21 into Eguation 2.20 and simpli- 
fying yields 

V p = SQRT[ (GM/Rp) ( He) ] (2.2 2) 

The required change in velocity can now be computed 
by subtracting V c from V p . 

V p - V c = SQRT[ (GU/Bp) (He)] - SQET[GM/E C ] (2.23) 

At perigee R c = R p ; V pc will represent the difference in 
vel ocit y. 

V pc = SQET (GM/£ C ) [ SQRT (He) - 1] (2.24) 

The square root of (1 + e) can be approximated by 
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1 * . 5e + . . - . Also note that SQRT (GM/R C ) = V c . Thus 

Equation 2.24 can be written as 

V pc = V c [ 1 * .5e - 1] (2.25) 



= V c e (2.26) 

2 

= v c [ (K+ 1-u/2tt) 2/3 - 1 ] (2.27) 

2[ K+l-u/2 tt ]2 / 3 

Equation 2.27 is a key formula and can be used to 
determine how much additional velocity is required to ensure 
a proper rendezvous transfer orbit. It has the flexibility 
to vary the number of circular orbits which is equivalent to 
being able to predetermine the rendezvous time. 

The use of this equation can be extended to 

deploying numerous, equally spaced satellites from a single 
large parent vehicle. The parent vehicle would be in the 
elliptical orbit carrying the satellites to be deployed. 

Each time the "rendezvous point" is achieved a satellite 
would be ejected. Two satellites would not De placed in the 

same position as the new rendezvous point is always u 

radians ahead of the current one. Because Equation 2.27 has 
u as a variable parameter, any angular positioning is 
possible and numerous satellites may be deployed limited 
only by the total weight tnat can be launched from Earth. 

3 • Orbital T ra n sfer Fuel Requireme nts 

The next parameter that must be calculated is how 
mucn fuel is required to inject the parent satellite from 
the circular orbit into the elliptical orbit. Comparison of 
Equations 2.7 and 2.22 show that V p has a larger magnitude 
tnan V c even though both are at a distance R p from the 
center of the Earth. Thus, a transfer from a circular to 
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elliptical orbit requires the satellite to accelerate in the 
direction of travel. Similarly, when the parent satellite 
ejects a payload, that payload must burn fuel and slow down 
to V c or it would continue on the parent's elliptical path. 
The amount of fuel required for either maneuver is critical 
since not only the fuel but also its container subtract from 
the weight that could be used for instrumentation. 

Rocket fuel energy is classified by its specific 
impulse (I sp ) and is a measure of the thrust developed as a 
function of fuel mass flow rate. The mass flow rate can be 
in terms of gallons per minute for a liquid fuel or pounds 
per second for a burning solid fuel. The accepted English 
units for I sp are expressed as seconds. This comes from 

X ^ p = pounds of thrust = seconds (2.28) 

pounds of fuel/second 

Typical values of I sp range from 275-325 seconds. Knowing 
the fuel I sp , the force of acceleration or thrust (F f ) can 
be found by 

F t = I sp (dW/dt) (2.29) 

W is the propellent weight flow rate. 

The rocket thrust is generated from burning fuel 
that exits the rocket at a velocity V x with units of feet 

pec second. The fuel I sp and exhaust velocity can De 

related by 

V x = gl sp (2.30) 

where g is the acceleration constant equal to 32.2 ft/sec 2 . 
The thrust can also be related to the exhaust velocity and 
fuel flow rate by 

F t = -(dm/dt) V x (2.31) 
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The minus sign is necessary because the exhaust is traveling 
in one direction and the resultant force is in the opposite 
direction. By Newton's laws, the rocket must react to the 
fuel thrust by 



F = a s a (2.32) 

d s is the mass of the rocket being accelerated and a is the 
resultant acceleration. The acceleration can be written in 
terms of the rocket velocity (V r ) as 

a = dV r /dt (2.33) 

Substituting Eguation 2.33 into Equation 2.32 gives 

F = M s (dV r /dt) (2.34) 

F t must equal F. Setting Equation 2.34 equal to Equation 
2.3 1 gives 

M s (dV r /dt) = - ( dm/ dt) V x (2.35) 

Simplifying and rearranging terms gives 

dV r = -V x (dm/M s ) (2. 36) 

The left side of this equation can now be integrated from 
the initial rocket velocity (Vj ) to the final rocket 
velocity (Vj ) . The right side can oe integrated from the 
initial rocket weight (W ; ) to the final rocket weight (Wf). 




Eguation 2.38 shows the evaluation result. 
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-V x ln (W f /W; ) 



(2. 38) 



V, -7, 



Dividing by -V x and using both sides as a power of the expo- 
nential e gives 



exp£-(V f -Vj )/V x ] = W f /W; (2.39) 

W, = Wj exp[ - (V f -V; ) /V x ] (2.40) 

The initial weight must be the sum of the final weight plus 
the weight of any fuel or propellent that was burned (W p ) to 
produce the change in velocity, or 

W f = Wj -W p (2.4 1) 

Substituting Equation 2.41 into Equation 2.40 gives 

Rj -W p = W; exp[- (Vf-Vj ) /V x ] 

Solving for W p gives 

R p = Wj [ 1 - exp- (V,-Vj ) /V x ] 

Replacing V x with gl sp from Equation 2.30 gives 

Wp = Rj r 1 - exp- (V, -Vj ) /(gl sp ) ] (2.44) 

The quantity Vf-Vj is the same change in velocity required 
by Equation 2.27 for orbital transfer. 

Vpc = Vf-Vj (2.45) 

Equation 2.44 then becomes 



(2. 42) 



(2.43) 
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w p = Wj[ 1 - ex p- (V pc /gl sp ) ] 



(2.46) 



Equation 2.27 can be used to determine the amount of 
velocity change that is needed for an orbital transfer and 
Equation 2.46 is used to determine the percentage of initial 
satellite weight that must be reserved for the fuel to 
accomplish the same transfer. 
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. SAT ELLITE DESIGN 



A. GENEBAL DESIGN 

The satellite will be placed in a circular low Earth 
orbit (LEO) of 322 km. An orbit close to the Earth will be 
used because of the anticipated low power level of the 
earthquake electromagnetic signal. The design of the satel- 
lite will be based on using the space shuttle for achieving 
the initial circular orbit. Because of the large cargo bay 
on the shuttle (60 ft X 15 ft) , a group of four satellites 
arranged in a circle and mounted on a central support bus 
will comprise the launch package. When deployed, the 

satellites will be egually spaced by 90 degrees. 

Primary electrical power will be provided by solar cells 
with batteries supplying power during the time the satellite 
is in the Earth’s shadow. Spin stabilization of the satel- 
lite body with hydrazine thrusters for station keeping will 
be used to maintain a de-spun antenna pointing at the Earth 
with attitude sensing to be done by Earth limb sensors. 
Each satellite will weigh approximately 500 pounds and have 
a cylindrically shaped main body. 

B. ANTENNA DESIGN 

The satellite will be designed to detect electromagnetic 
energy in the 30 MHz region. This frequency is used because 
it is designated for astronomical observations and manmade 
interference should be at a minimum. Additionally 30 MHz is 
well above the cutoff frequency (F p ) for propagation through 
the ionosphere given by 

F p = 9SQRI ( N e ) (3.1) 

where N e is the ionospheric plasma electron density with 
units of electrons/meter 3 . A nominal value for F p is 
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10 MHz. Any frequency higher than F p will successfully 
penetrate into outer space with lower frequencies being 
reflected back to Earth. 

A Yagi antenna will be used as the main sensor because 
of the size and space constraints imposed by a satellite 
design. Reference 8 will be used as a basis for the design 
of the antenna. Table I [Ref. 8:p. 7] is used once the 

initial design criteria is established. This initial data 
consists of the frequency or operation (30 MHz) , antenna 
gain desired (9.2 dB ) , d/A ratio where d is the diameter of 
the parasitic elements (2 cm) , and the D/a ratio where D is 
the diameter of the support boom (4 cm) . The wavelength A 
is given by 

A = c (3.2) 

f 

The speed of light (c) is a constant egual to 
3 X 10 8 meters/second and f is the operating frequency. 
Thus a frequency of 30 MHz gives a wavelength of 10 meters. 
Following the steps in Reference 8, Figure 3.1 [Ref. 8:p. 
20] gives the length of the first and third director as 
4.44 meters, 4.5 meters for the length of the second 

director and 4.9 meters for the length of the reflector. 
The element spacing is 2 meters. Deployment of the antenna 
in outer space can be achieved by using a bistem which is a 
deployable boom rolled up on a drum in its flattened config- 
uration. Because of its small size when stowed, total 
storage space required will be approximately one-half a 
cubic foot. Since the antenna will be operating in a 

weightless environment, strong structural members are not 
required. The total weight of the antenna would not exceed 
30 pounds since the material of construction could be a 
light-weight aluminum based alloy. The gain of 9.2 dB is 
referenced to a dipole antenna. For reference to an 
isotropic radiator, 2.16 dB of gain must be added. 
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TABLE I 

YAGI ELEMENT LENGTHS 
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C. FUEL BEQUIBEfl EHT S 
1 . Sat ellite Fuel 

The on-station weight of the satellite is 500 pounds 
and this weight must be decelerated from the parent* s ellip- 
tical velocity. Equation 2.27 will be used to determine the 
necessary velocity change and Eguation 2.7 will be used to 
find V c . G!1 evaluates as a constant equal to 

3.99 X 10 14 m 3 /sec 2 and r is the sum of the radius of the 
Earth {6370 Km) and the satellite altitude above the 
Earth (322 Km) . 



V c = SQET (3 . 99 X 1 0 14 /669 2 X 1Q3) 



(3-3) 



= 7,722 m/sec 



(3.4) 



= 25,335 ft /sec (3.5) 

Equation 2.27 can now be evaluated if K is set equal to one 
and u is set equal to ir/2 radians. 

V pc = 25 , 33 5[ (2-. 25) 2 /3 - 1] (3.6) 

2[ (2-. 25) 2/3 - i ] 

= 3,945 ft/sec (3.7) 

Equation 2.46 can be used to find the fuel fraction given by 
Wp/W ; = 1 - exp-[ 3945/(32) (300) ] (3.8) 



= .337 



(3.9) 
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Figure 3.1 lagi Antenna Design Curves 
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ELCHtNT OIMETEA TO WAVELENGTH PATIO, i/X 



The initial weight is composed of 


the cn- orbit 


satellite weight of 500 pounds (W s ) , W p and 
the case (W c ) that contains the fuel. That is 


the weight of 


Wj = Hs*Wp + Kc 


(3. 10) 



For an efficiently designed rocket motor, the majority of 



the weight is due to fuel with the remaining 
the case. This ratio (f) is given by 


weight due to 


f = w p 

«p*W c 


(3. 11) 



A typical value of f is 0.93. Using this definition of f. 
Equation 3.10 can be rewritten as 

H; = W s +W p /f (3.12) 

Substituting H p = . 3 3 7 W j from Equation 3.9 and 0.93 for f 
into Eguation 3.12 gives 

Hi = W s + .337W; /. 93 (3. 13) 



Rearranging terms to solve for Wj and 

500 pcunds for W s gives 


subs titutin g 


Wj = 500 

1 - .337/. 93 


(3. 14) 


= 784 lb 

This can be used to find W p from Equation 3.9 


(3. 15) 


W p = (784) (.337) 


(3. 1b) 
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264 lb 



(3. 17) 



W p = 



The weight of the case 
Equation 3.11 



can be found by rearranging 



H c = Wp(l-f) 



{3. 18) 



= 264 (1-.93) 
. 93 



(3. 19) 



= 20 lb 



(3. 20) 



2 . Parent Sa tellite Fu el 

The amount of fuel required by the parent to inject 
itself into the elliptical orbit can now be calculated in a 
similar manner. The central support bus weight will be set 
at 500 pounds. This relatively large weight is required 
because the bus must be strong enough to rigidly hold the 
four satellites during the sh uttle- launch. Thus, enough 
fuel must be allowed to accelerate the bus and three satel- 
lites into an elliptical orbit since one satellite is 
deployed prior to the parent's acceleration. The net weight 
to be boosted is then 2,852 pounds. Similar to Equation 
3.14, the initial total weight is 



W; = 2852 



(3.21) 



1 - .337/. 93 



= 4,473 lb 



(3. 22) 



The fuel weight and case weight for the parent can also be 
found in the same manner as for the individual satellite. 



Wp = (4473) (.337) 



(3. 23) 
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1,507 lb 



(3.24) 



H 



P 



W c = 1507 (1-. 93) 
. 93 



(3.25) 



= 114 lb 



(3. 26) 



D. ELECTRICAL POWER 

The primary source of electrical power will be provided 
by converting light energy into electrical energy by use of 
solar cells. Because the available power from an individual 
cell is small (on the order of 0.06 watts), a large number 
of cells must be used that collects solar energy over a 
large incident area. In order to provide the desired total 
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Figure 3.2 Satellite Power System 



vary as a function of load, temperature and other factors. 
Bus voltage limiters prevent the voltage swings from 
exceeding spacecraft design ranges by acting as variable 
shunt loads during a rising voltage condition. The battery 
discharge regulator turns on when a pre-determined bus low 
voltage is sensed and provides a regulated output. The 
battery terminal voltage is higher than the required bus 
voltage and in order to charge the battery a still higher 
level of voltage is needed. To provide this higher voltage 
level, the additional solar panel is series connected to tne 
main panel. The size of the battery charging panel is 
designed so that the charge control functions provided in 
the controller need only be an on-off and rate change 
switch. The voltage and current sensors provide means to 
monitor the operational state of the power system and for 
load managing during the orbit. 

1 . Solar Pan el Design 

As shown in Figure 3.3 [Ref. 9:p. 1-7.11] individual 
solar cells are arranged in a series-parallel array to 
provide the required power. The three strings are in 

parallel to minimize power loss should a single cell fail. 
These arrays are mounted on the satellite with the series 
string length parallel to the satellite spin axis. All the 
arrays are then connected in parallel to produce the desired 
power. The surface of the solar cell is considered to be a 
Lambertian surface and the power produced by the cell 
decreases by a cosine function as the sun's incident angle 
changes from normal as shown in Figure 3.4 [Ref. 9:p. 
1-7. 11]. 

Figure 3.5 [Ref. 9 : p. 2-4.17] shows how the sun's 
incident angle varies during one year. When the solar cells 
are in the sunlight, only those positioned in the center are 
illuminated by normal rays. Going left or right from center 
causes the ceil output to drop by the cosine function. The 
net result is an output from the half-cylinder that is equal 
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Figure 3.3 Solar Cell Array 

to the output from a flat plate whose dimensions are the 
diameter and height of the cylinder. Additionally, the sun 
line varying 23.5 degrees as shown in Figure 3.5 causes a 
further decrease in solar panel output power. 

Solar cell performance also degrades with time from 
exposure to electron and proton fields trapped in the 
Earth's magnetic field and solar flares. The extent of 
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Figure 3.4 Cylindrical Panel Installation 

degradation is shown is Figure 3-6 [Ref. 9: p. 1-7.13]- 
Damage is a result of charged particles creating defects in 
the n and p layers of the cell. The net result is shown in 
Figure 3.7. For any set of four sequential peaks, the maxi- 
mums correspond to the autumnal and vernal equinox and the 
minimums correspond to summer and winter solstice when sun 
incidence is offset by 23.5 degrees. 



43 




Figure 3.5 Solar Incident Angle 

The satellite that is being proposed will have a 
diameter of 5.7 feet and a height of 10 feet. This gives a 

projected area of 57 square feet. At beginning-of-lif e 
(BOL) solar cells produce approximately 11 watts of power 
per sguare foot of area. Thus the 57 square feet will yield 



44 



r 



1 




Figure 3.6 
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Figure 3.7 Panel Output VS Time 



electrical power budget must be based on this EOL value. 
The p over- to-mass ratio at EOL for solar cells is approxi- 
mately 20 watts/kg for normal incidence [Ref. 10:p- 73]. 

This ratio must be decreased because the satellite does not 
have flat plate collectors. The amount of decrease will be 
the ratio of the area of a flat plate collector to the area 
of the illuminated half cylinder. The area of an equivalent 
flat collector is 57 square feet. The surface area of the 
half cylinder (A c ) is 

A c = (.5) (277-rh) (3.28) 



= (*) (5.7/2) (10) 



(3. 29) 



= 89.5 ft* 



(3.30) 



The power- to- mass ratio of the cylinder (P c ) is then 



P c = (20) (57/89.5) 



(3.31) 



= 12.73 watts/kg (3.32) 

The mass of the solar cells (M c ) is then 

M c = 200/12.73 (3. 33) 

= 15. 7 kg (3. 34) 



This is the mass of only those cells that are receiving 



sunlight. The total mass is twice as much. Converting kg 

to pounds will give a total weight of the solar cells as 
69 pounds. 

2 - Batter y R equirem ents 

Batteries that are used in a weightless environment 
can not he of the "flooded" electrolyte type such as lead- 
acid cells. Rather, the electrolyte must be absorbed in the 
separator material to assure even distribution. The two 
primary types of batteries that meet this condition are the 
sealed silver-zinc and the nickel-cadmium. 

Table II [Ref. 9:p. 1-7.17] shows the major charac- 

teristics of concern for both battery types. Nickel-cadmium 
batteries will ne used because of their longer life expec- 
tancy and tolerance to overcharge allowing a simpler charge 
control circuit. The major disadvantage will be an increase 
in weight due to the lower energy density. 

The satellite will be spending 50 percent of its 
time in the Earth’s shadow during which tne battery must 
supply the electrical power. During the time sunlight is 
available, the battery must be recharged. Since equal 

amounts of time are spent in light and shade, half of the 
solar panels must be dedicated to battery charging and the 
remaining half to satellite bus power. The time the satel- 
lite is in darkness (T d ) can be found from Equation 2. 10 to 
be 



= . 5T C 



(3. 35) 



= (tt) (6692 X 10 3) 3/2/sqrt (3.99 X 10**) 



(3. 36) 



T d = 4 5. 4 min 



(3.37) 
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TABLE II 

BATTEBY CHARACTERISTICS 
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200 watts of power 
required from the 



The satellite system will be drawing 
during this time and the watt-hours 
battery (P b ) will be 

P b = (200) (45.4/60) (3.38) 



= 151.3 watt-hours (3.39) 

Figure 3.8 [Ref. 9:p. 1-7.21] shows tnat as the 

depth of discharge (DOD) decreases the lifetime of the 
battery increases. In order to have a long lifetime, a high 
capacity battery system is desirable but the penalty is 
additional weight. From Reference 10, the battery system on 
the Intelsat V can supply 34 amp-hours at 27 volts and has a 
mass cf 32 kq [Ref. 10:p. 73]. This gives a battery energy 
density of 13 watts/pound which is slighty better than shown 
in Table II. To provide a low DOD but not require excessive 
weight, the battery system for the proposed satellite will 
be capable of supplying 1200 watt-hours. This gives a DOD 
of 



DOD = 15 1.3/1200 



(3.40) 



= 12.o percent (3.41) 

From Figure 3.3 this gives approximately 23,000 cycles for 
the battery. The time for one cycle is 90.8 minutes. There 
are 5,789 cycles completed in one year. Thus the battery 
has a lifetime of approximately four years and will weigh 
about 92 pounds. 
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Figure 3-8 Nickel-Cad aium Battery Cycle Life 
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CELL TEMPERATURE, °F 



E 



SYSTEM WEIGHT AND POWER BUDGET 



The satellite has an on 
Subtracting 30 pounds for the 
solar cells and the battery 
309 pounds for the remaining 
control systems. Using the 
Reference 10 gives the results 



orbit weight of 500 pounds, 
antenna, 69 pounds for the 
weight of 92 pounds leaves 
electronics, structure and 
weight and power ratios in 
shown in Table III. 



TABLE III 

WEIGHT AND POWER BUDGET 



System 


Height (lb) 


Power ( watt; 


Antenna 


30 


— 


Solar 

Cells 


69 




Battery 


92 


— 


Telemetry , 
Command 


16 


7.9 


Attitude 

Control 


43 


9.6 


Propulsion 


24 


1 . 0 


Structure 


86 


— 


Electrical 

Integration 


25 


20.3 


Signa 1 
Processin g 


1 15 


161.2 



Of the 161.2 watts available for signal processing, 
approximately 30 percent would be available for use to 
transmit data back to central Earth collection stations. 
The signal strength at the surface of the Earth can be 
calculated by 

p r = p t (G t G r ) A 2 (3.42) 

(4tt) 2 r 2 
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